February 3, 2008 8:15 



WSPC - Proceedings Trim Size: llin x 8. Sin osullivan 



1 



Parsec-scale Investigation of the Magnetic Field Structure of Several AGN Jets 
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Multi-frequency (4.6, 5, 5.5, 8, 8.8, 13, 15, 22 &: 43 GHz) polarization observations of 6 "blazars" were obtained on 
the American Very Long Baseline Array (VLB A) over a 24-hr period on 2 July 2006. Observing at several frequencies, 
separated by short and long intervals, enabled reliable determination of the distribution of Faraday Rotation on a 
range of scales. In all cases the magnitude of the RM increases in the higher frequency observations, implying that the 
electron density and/or magnetic field strength is increasing as we get closer to the central engine. After correcting 
for Faraday rotation, the polarization orientation in the jet is either parallel or perpendicular to the jet direction. A 
transverse Rotation Measure (RM) gradient was detected in the jet of 0954+658, providing evidence for the presence 
of a helical magnetic field surrounding the jet. For three of the sources (0954+658, 1418+546, 2200+420), the sign 
of the RM in the core region changes in different frequency-intervals, indicating that the line-of-sight component of 
the magnetic field is changing with distance from the base of the jet. We suggest an explanation for this in terms of 
bends in a relativistic jet surrounded by a helical magnetic field; where there is no clear evidence for pc-scale bends, 
the same effect can be explained by an accelerating/decelerating jet. 
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1. Introduction 

AGN jets emit synchrotron radiation that can be de- 
tected at ah radio frequencies. This radiation is often 
highly hnearly polarized, which provides important 
information on the degree of order and orientation 
of the magnetic field in these jets. The type of AGN 
studied in this experiment are known as "blazars", 
which have jets pointed close to our line of sight 
(LoS) and often exhibit strong variability in total 
fiux and linear polarization over a broad range of 
frequencies from 7-ray to radio. 

We observed 6 sources with the American Very 
Long Baseline Array (VLBA) at 8 frequencies from 
4.6 GHz to 43 GHz over a 24-hr period on 2 July 
2006. This radio interferometer provides milliarc- 
second resolution which corresponds to the parsec 
scale structure of these jets. Even though these AGN 
are intrinsically two-sided, pc-scale observations typ- 
ically show a one sided "core-jet" structure due to 
Doppler boosting of the radiation from the relativis- 
tic jet pointed towards us. The jet moving away from 
us is usually too faint to be detected. 

The radiation detected from the jet is gener- 
ally optically thin while the core region usually dis- 
plays a fiat spectrum attributed to synchrotron self- 
absorption. In the optically thin regime the polar- 
ization orientation is perpendicular to the magnetic 
field, while in the optically thick/fiat spectrum region 



the polarization is parallel to the magnetic field. 

Recent MHD simulations have provided an al- 
most complete explanation of how these jets are 
launched, accelerated and collimated close to the 
black hole, see [1] and references therein. The global 
magnetic field structure expected on these scales is 
helical. It is possible that this structure changes when 
the flow gets disrupted by shocks after a few hundred 
Schwarzschild radii, but observational evidence of he- 
lical fields on scales larger than this [2-4] suggests 
that remnants of the earlier magnetic field structure 
remain or that a current driven helical kink instabil- 
ity is generated [5, 6]. 

2. Faraday Rotation 

A rotation of the plane of polarization of an electro- 
magnetic wave occurs when it propagates through a 
region with free electrons and a magnetic field, com- 
monly referred to as Faraday Rotation. The effect 
is manifest as a linear dependence of the polariza- 
tion angle (x) on the wavelength (A) squared and is 
described by the formula 



Xobs = Xo + 



27rm^c' 



(1) 



where the Rotation Measure (RM), in parentheses, 
is proportional to the integral of the electron den- 
sity and the dot product of the magnetic field along 
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the LoS with the path length dl. Hence, a posi- 
tive/negative RM tells you that the LoS magnetic 
field is pointing towards/away from the observer. 

For a pure e+/e~ plasma one would expect 
no Faraday rotation because of cancellation due to 
the dependence. The dependence on the inverse 
square of the mass also shows that the effect is much 
greater for thermal electrons than for protons or rel- 
ativistic electrons. Therefore, multi-frequency obser- 
vations enable us to see gas that would otherwise 
be invisible to us and, combined with the corrected 
polarization orientation, we are provided with a 3-D 
view of the magnetic field structure. 



1156+295 



+136rad/m^ 




Fig. 1. Plots of X vs. for the whole observed frequency 
range of the core regions of 1156+295 (right) and 0954+658 
(left) clearly showing a transition in both cases where physi- 
cally different regions are being probed 



However, if the rotating medium is mixed in with 
the emitting plasma then internal Faraday rotation 
occurs and the above formula does not generally ap- 
ply. Theoretical models [7, 8] of internal FR in a 
spherical or cylindrical region have shown that, if a 
rotation of greater than 45 degrees is observed, the 
Faraday rotation must be external. Many observa- 
tions show rotations greater than 45 degrees, [3, 9, 
10] and our current observations, indicating that the 
rotating gas must be segregated from the emitting 
region. Furthermore, observations of short term RM 
variability [11] means that the Faraday rotating gas 
must be very close to the source, most likely in a 
boundary layer or sheath surrounding the jet. 



3. Observations 

The calibration and imaging was performed with the 
radio astronomy software package AIPS. The inte- 
grated foreground (mainly Galactic) RMs were sub- 
tracted using the multi-frequency VLA observations 
of [12]. This isolates the RM distribution in the im- 
mediate vicinity of the AGN and is very important 
for source rest frame calculations, since the observed 
RM scales with (1 + z)^. 

Our current observations are ideally suited for 
Faraday rotation analysis, with both long and short 
spaced frequency intervals, which significantly re- 
duces the ztnTT ambiguity in the observed polariza- 
tion angles. This enables us to construct RM maps 
on a range of different scales and Faraday depths 
providing a wealth of information on the RM dis- 
tribution and sign in the core region and out along 
the jet. Matched-resolution images were constructed 
for different frequency intervals where clear transi- 
tions in RM occurred, corresponding to sampling of 
physically different regions (eg. Fig. 1). 




MilMARC SEC 

Fig. 2. The left image of 2200+420 clearly shows how the po- 
larization sticks are aligned with the jet direction and, impor- 
tantly, remain aligned even as the jet bends. The image on the 
right of 1418+546 is the one source that displays polarization 
perpendicular to the jet direction. The contours correspond to 
total radio flux while the sticks correspond to the polarization 
orientation. 

4. Results 

Since Faraday rotation has a very significant effect 
on the observed polarization angle, Xobs^ the rotation 
must be removed in order to obtain the intrinsic po- 
larization orientation, xo^ and reach any definite con- 
clusions about the magnetic field structure. From our 
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multi-frequency observations we find that the polar- 
ization in the optically thin jet is either aligned with 
(0954+658, 1156+295, 2007+777, 2200+420) or per- 
pendicular to (1418+546) the observed jet direction 
(eg. Fig. 2). ^ The aligned polarization cases can be 
attributed to the toroidal component of a relatively 
tightly wound helical field, whereas the poloidal com- 
ponent, which dominates in a loosely wound helical 
field, can account for the one perpendicular polariza- 
tion observation [13]. 



Table 1. Note the higher magnitude RM in the high u 
range for all sources, as well as the RM sign change in 
half the sources. All RM values in rad/m? 



Blazar 


z 


Core RM 


Core RM 






(Low u range) 


(High u range) 


0954+658 


0.368 


+33 ± 14 


-1591 + 265 


1156+295 


0.729 


+136 + 4 


+1647 ± 159 


1418+546 


0.152 


+75 + 8 


-501+48 


1749+096 


0.320 


-33 ± 24 


-500 ± 347 


2007+777 


0.342 


+638 ± 39 


+1946 ± 140 


2200+420 


0.069 


+240 ± 90 


-1008 + 43 



In all cases, the RM has a larger magnitude in 
the core region than in the jet for a particular fre- 
quency interval. The magnitude of the core RM also 
increases in the higher frequency intervals for all 
sources, simply implying that the electron density 
and/or magnetic field strength is increasing as we 
sample regions closer to the central engine. This can 
be clearly seen in Table 1 where a summary of the 
core RM results is presented. The redshift, z, is also 
included. 

Head-on (0 < l/F) Tail-on (0 > \IT) 

^ "^ ^^^^^^d^ ^ Direction | ""^»>...^^9 Direction 

\ Viewed at < 1/r / ^^^^5^ 

\ in observer's frame / Viewed at > 1/r 

^ ^ in observer' s frame 




Fig. 3. For a "head-on" {0 < l/F) or "tail-on" {$ > l/F) 
view of a hehcal magnetic field, different LoS components will 
dominate (shown by solid circles surrounding the dot or cross 
which indicate the direction of the LoS component) producing 
RMs with different signs in an unresolved jet. 



^No appreciable jet polarization was detected in 1749+096 



A surprising result is that, for three of the 
sources, the sign of the RM changes in different fre- 
quency intervals, telling us that the LoS component 
of the magnetic field is changing with distance from 
the base of the jet. In the case of 2200+420, the 
dominant jet magnetic field is transverse through- 
out the jet, even as it bends (see Fig. 2), consistent 
with the presence of a global helical magnetic field 
structure. Bends in a jet, surrounded by a helical 
magnetic field, on scales smaller than is probed by 
the lowest frequency resolution will produce differ- 
ent dominant LoS magnetic field components due to 
the relativistic motion of the jet towards us, see Fig. 
3. This effect provides a natural explanation of the 
observed RM sign changes. 




Fig. 4. RM maps of 1418+546 in different frequency inter- 
vals showing how the sign of the RM changes in the core. Grey 
scale denotes the RM in rad/m? and contours are of total ra- 
dio flux. The insets are plots of x vs. for the corresponding 
regions. 

If there is no clear evidence for bends in the in- 
ner jet, as is the case for 0954+658 and 1418+546 
(eg. Fig. 4), then the same effect can occur for an 
accelerating/decelerating jet surrounded by a heli- 
cal magnetic field; although it remains possible that 
bends are present on scales smaller than our resolu- 
tion. 

Further analysis of 0954+658 reveals the pres- 
ence of a RM gradient across the jet (Fig. 5), which 
is a strong signature for the presence of a helical mag- 
netic field geometry. For a side-on view of a helical 
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magnetic field, one would expect a systematic gradi- 
ent across the jet with the LoS magnetic field chang- 
ing sign. However, as was shown in [2], a same-sign 
gradient (as detected in 0954+658) can be explained 
by having the viewing angle {0) in the jet rest frame 
smaller than the pitch angle {ip) of the helix. 

Another excellent diagnostic of the magnetic 
field structure are transverse profiles of the jet. For a 
helical magnetic field with < t/j^ one would expect 
the peak in the polarization profile to be offset from 
the peak of the total flux profile, and furthermore, 
the magnitude of the RM gradient should increase 
on the side opposite to the total flux peak offset [14]. 
Transverse profiles of the jet of 0954+658 display this 
pattern, see Fig. 5. 



Total Flux|[ Polarized Flux]| Rotation Measur^ 




MilliARC SEC 



Fig. 5. The grey scale image displays the RM gradient in 
0954+658 with accompanying scaled transverse profiles of to- 
tal fiux (solid line), polarized fiux (dotted line) and RM inten- 
sity (grey line) across the slice indicated. 

5. Conclusions 

Several phenomenon, for example, the polarization 
structure [13], RM gradients [2] and circular polar- 



ization generation [15], can be understood by consid- 
ering a helical magnetic field geometry for pc-scale 
AGN jets. In this paper, we consider how regions 
with different RM signs can also be explained within 
a helical magnetic field model, as places where the 
jet is observed at angles greater than or less than 
1/r due to bends in the jet or due to an acceler- 
ating/decelerating straight jet. (A longitudinal jet 
magnetic field with a change in the angle to the 
LoS could also cause a RM sign change, but this 
does not correspond to the observed magnetic field 
in 0954+658 or 2200+420.) It's important to note 
that VLBI resolution is usually not sufficient to com- 
pletely resolve the true optically thick core, there- 
fore, the VLBI "core" consists of emission from the 
true core and some of the optically thin inner jet. So 
if bends occur on scales smaller than the observed 
VLBI "core", core RMs with different signs can be 
derived from observations at different frequencies (ie. 
probing different scales of the inner jet). 

In our future work, we will attempt to recon- 
struct the 3-D path of the jet through space us- 
ing the combined information from the observed dis- 
tributions of the total intensity, linear polarization, 
spectral index and rotation measure. 
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